Magnetoplasmons in layered graphene structures 
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We calculate the dispersion equations for magnetoplasmons in a single layer, a pair of parallel 
layers, a graphite bilayer and a superlattice of graphene layers in a perpendicular magnetic field. We 
demonstrate the feasibility of a drift-induced instability of magnetoplasmons The magnetoplasmon 
instability in a superlattice is enhanced compared to a single graphene layer. The energies of the 
unstable magnetoplasmons could be in the terahertz (THz) part of the electromagnetic spectrum. 
The enhanced instability makes superlattice graphene a potential source of THz radiation. 

PACS numbers: 71.35.Ji, 71.35.Lk, 71.35.-y 



I. INTRODUCTION 

Recent advances in fabrication techniques have made it 
possible to produce graphene, which is a two-dimensional 
(2D) honeycomb lattice of carbon atoms forming the ba- 
sic planar structure in graphite [l:]. Graphene has stim- 
ulated considerable theoretical interest as a semi-metal 
whose electron effective mass may be described by an 
unusual massless Dirac-fermion band structure. Several 
novel many-body effects in graphene have been inves- 
tigated 0, Q. The theory of Weiss oscillations in the 
magnetoplasmon spectrum of Dirac electrons in graphene 
has been developed in Ref. [ij . The magnetoplasmon ex- 
citations in graphene for filling factors v < 6 has been 
calculated in Ref. ||. In recent experiments, the inte- 
ger quantum Hall effect (IQHE) has been reported @. 
Quantum Hall ferromagnetism in graphene has been in- 
vestigated from a theoretical point of view [7| . Graphene 
has a number of interesting properties as a result of its 
unusual band structure which is linear near two inequiv- 
alent points (K and K') in the Brillouin zone. In the 
presence of a magnetic field, the graphene structure shifts 
both the Shubnikovde Haas oscillations 0] as well as the 
step pattern of the IQHE [9(. Both these effects have 
recently been reported experimentally [6| . The spectrum 
of plasmon excitations in a single graphene layer embed- 
ded in a material with effective dielectric constant £& in 
the absence of an external magnetic field, has been calcu- 
lated in In this paper, we show that features, such 
as charge density oscillations, arise when a magnetic field 
is applied. 

This paper is organized as follows. In Sec. [Til we ana- 
lyze the magnetoplasmon spectrum for a single graphene 
layer. The collective charge density excitations in a 
strong magnetic field for a graphite bilayer and a bilayer 
graphene are calculated in Sees. IIIII and IIV[ respectively. 
The enhancement of the magnetoplasmon instability in 
an infinite periodic graphene superlattice is investigated 
in Sec. |V] In these calculations, we assume that there is 
no tunneling between the graphene layers forming the su- 



perlattice. The results of our numerical calculations are 
presented in for each structure investigated. A brief dis- 
cussion of plasmon instabilities in graphene is presented 
in Sec. ED 



II. A SINGLE GRAPHENE LAYER 

Let us first consider electrons in a single graphene layer 
in the xy-plane in a perpendicular magnetic field B par- 
allel to the positive z axis. Here, we neglect the Zeeman 
splitting and assume valley energy degeneracy, describ- 
ing the eigenstates by two pseudo spins d, 03 • We have 
an effective 2x2 matrix Hamiltonian H/gs whose diag- 
onal elements are zero and whose off-diagonal elements 
are n x ± iiTy where tt = —iKV + eA, — e is the electron 
charge, A is the vector potential, vp — \/3at/(2H) is the 
Fermi velocity with a = 2.566A denoting the lattice con- 
stant, and t w 2.71 eV is the overlap integral between 
nearest- neighbor carbon atoms Q. 

Choosing A = (0, Bx, 0), the eigenfunctions of 7J( ) 
are labeled by a = {k y , n, s(n)}, where n = 0, 1, 2, • • • 
is the Landau level index, k y is the electron wave vec- 
tor in the y-direction, and s(n), which is defined by 
s(n) = for n — and s(n) = ±1 for n > 0, 
labels the conduction (+1) and valence (—1 and 0) 
band, respectively. The eigenfunction V'a(r) is given 
by a spinor ip a (x,y) with components given by Q 

!,(2) _ 



i(x + lj I ky)/y/L y and ipa' 
L y . Here, Ijj = y/H/eB, 



= C„e lfc ^s(n)i"- 1 $„_ 

C n e ik yyi n $ n (x + l%ky)/ 
and L y is a normalization length. We have C n = 1 
for n = 0, C n = l/\/2 for n > and <J> n (x) = 

(2"n!V^r 1/2 e^ x ' l ^ 2 l 2 H n {x/l H ), where H n (x) is 
a Hermite polynomial. The eigenenergies are given 
by e a = s(n)e n = s(n)(hvF /Ih)V^i, for which successive 
levels are not equally separated. 

The dynamic dielectric function in RPA [13| is given by 
e(q,uj) = 1 — V c (q)H(q, u>), where q is the in-plane wave 
vector, V c (q) = 2ite 2 / (e s q) is the 2D Coulomb interaction 



2 



and the 2D polarization function is 
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where f s i n ) n is the Fermi-Dirac function, i^ ss / (n, n') 
arises from the overlap of eigenstates and is given by 
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The magnetoplasmon dispersion relation for a single 
graphene layer was obtained by seeking the solutions of 
e(q,w) = 0. The highest valence band is full and all 
others empty at T = K. Transitions to the lowest five 
Landau levels in the conduction and valence bands were 
the only single-particle excitations included in our cal- 
culations. Fig. Q] (a) is the solution of the dispersion 
equation for a single layer of graphene when the imagi- 
nary part of the plasmon frequency is zero. In this case, 
the plasmons are self-sustaining oscillations except when 
they enter the particle-hole mode region where they un- 
dergo a loss due to Landau damping. In Fig. [1] (b), we 
plot the solutions of the dispersion when the frequency 
is complex for which the real part that is linear in q and 
exists only where the magnetoplasmon in Fig. Q] (a) has 
negative group velocity. The real and imaginary parts 
of the frequency are denoted by u>r and oji > 0, respec- 
tively. The loss corresponds to finite imaginary part of 
frequency. This instability after excitation is due to a 
transfer of energy back from a magnetoplasmon to an 
electric current which excites it, thereby making this col- 
lective mode unstable (see Fig [T] (b)). Thus, we have 
a non-zero imaginary part of the frequency for a single 
graphene layer in a magnetic field. The non-zero imagi- 
nary part for collective excitations for a 2D electron gas 
(2DEG) in semiconductors has been established, out only 
for several layers of semiconductor [I?], [H, [H, H(| • See 
also [11 [11 (Mill.) 

The negative group velocity for qln > 1 is caused by 
the magnetic field [21[. We used the same parameters 
are employed in calculating Fig. 1, but we summed over a 
larger number of Landau levels in the conduction and va- 
lence bands. Qualitatively, the results are the same. The 
main differences are that the number of single-particle 
excitation lines which are allowed is increased and the 
frequency of the highest mode which is, of course, af- 
fected by the number of Landau levels included in the 
sum. However, the lower branches of collective modes do 
not change significantly. Therefore, including in the cal- 
culations the five lowest Landau levels in the conduction 
and valence bands is justified. 




FIG. 1: Magnetoplasmon excitation energy as a function of 
wave vector, in units of , in a single graphene layer, (a) 
Real frequency solution, (b) The real and imaginary parts of 
the frequency satisfying the dispersion equation. 



III. A GRAPHITE BILAYER 

While the electron effective mass in a graphene single 
layer is zero, a graphite bilayer consisting of a pair of 
parallel graphite layers with inter-planar separation c/2 
implies finite electron mass. The electron spectrum a 
graphite bilayer in a magnetic field is very much different 
from the case of a single graphene layer. This is caused by 
interlayer hopping. Here, c/a = 2.802 with a = 2.566A 
denoting the lattice constant [22j]. The nearest-neighbor 
tight-binding approximation yields a gapless state with 
parabolic bands touching at the K and K' points in- 
stead of conical bands [23L l2~ij . A graphite bilayer can be 
treated as a gapless semiconductor. The eigenfunction 
V'a(r) of an electron in a graphite bilayer in a perpendicu- 
lar magnetic field is given for low-lying energy excitations 
by0 



C n } „ik„y ( ®n(x + l 2 H ky) 

s(n)Q n $ n - 2 {x + Ijiky) 



.(3) 



where a — {k yi n, s(n)}, Cn^ = 1 when n = or n = 1 
and Cn^ = l/\/2 when n > 2. Also, Q n — when 
n = or n = 1 and Q n — 1 when n > 2 and <fr n {%) is 
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defined above. The corresponding eigenenergy is e« = 
s(n)e^ = s{n)hio c yJn{n — 1), where lu c = eB/m with 
to = 7i/(2v 2 ), 71 = 0.39 eV and v = 8 x 10 5 m/s [H[ 
(compare to the electron spectrum in magnetic field in a 
single graphene layer presented above). 

Following the procedure described above, it can be 
seen that for the RPA dielectric function (q, u) for 
a graphite bilayer, we must replace the polarization by 
Il( b \q,uj) instead of H(q,ui). This is obtained from Eq. 



instead of 



(P) by means of the eigenspectrum e„ s / n ) 
£n,s(n)- I n this case, the form factor F ss i (n, n' , q) is re- 
placed by 



5 
4 

- 3 

LL 
> 

I 2 
1 




~i 1 1 r 



_i i i i_ 



,(6) 



1(n,n',q)=A(c^)' 
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dx exjp[iq x 3 



+ a(n)s'(n r )Q n Qn> / dx exp[iq x x]® n - 2 (x) 

J — oo 

x ^ n ^ 2 (x + l 2 H q y )\ 2 J . (4) 

The eigenfunction of bilayer Bernal graphene presented 
in Eq. ([3]) was obtained in Ref. [24| by taking into ac- 
count the interlayer Coulomb interactions whose dom- 
inant contributions are included. These are due to 
nearest-neighbor intralayer hopping (see Fig. 1 in Ref. 
|24| ) . The interlayer Coulomb interactions are included in 
Eq. (|4]) since they are contained in the wavefunctions 
of Eq. ([3]), entering into the overlap integral Eq. ((4|. 
In Fig. [21 we present the dispersion relation. The four 
straight lines correspond to the single-electron transi- 
tions between different Landau levels. Three curved lines 
are the undamped magnetoplasmon excitations. For a 
range of wave vectors the group velocity is negative due 
to the magnetic field, a result similar to those in a single 
graphene layer. The transfer of energy between collective 
excitations and electrons appears only when the charged 
particle velocity has the same value as the phase velocity 
of the collective mode. 
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FIG. 3: Magnetoplasmon dispersion relation for bilayer 
graphene at T — K with separation D = Ih- Only the high- 
est Landau level in the valence is occupied and completely 
full. 



qD 3> 1, Eq. §5§ reduces to the dispersion equation 
for magnetoplasmons in a single graphene layer. If 
Rn(q,u) = n 22 (<7,w) = U(q,w), and e x — s 2 — e b , then 
we get from Eq. ([5]): 



IV. A BILAYER GRAPHENE 

For bilayer graphene with layer separation, D and no 
interlayer hoping, we have the dispersion equation [ToL 
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smh 2 {qD) ( 2V c {q)nn(q,cj) - — - coth(qD) 
V e b 

x ^0(9^22(9,0;) - |J - coth(qD)^ = 1 , (5) 

where 11^(9, oj) is the polarization function of the charge 
carriers on the first j = 1 or the second j = 2 graphene 
layer defined by Eq. ((TJ. We observe that in the limit 



The dispersion relation for magnetoplasmons in bilayer 
graphene is presented in Fig. [3] These results show that 
each originally degenerate magnetoplasmon mode in each 
layer of an isolated single graphene layer is shifted from 
their value by the interlayer Coulomb interaction. For a 
range of wave vectors the group velocity is negative due 
to the magnetic field, analogous to the magnetoplasmon 
modes in a single graphene layer. A region of instability 
also exists for bilayer graphene. 
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FIG. 4: The imaginary part of the magnetoplasmon energy in 
a single graphene layer compared to the results for an infinite 
superlattice of graphene layers at v = 1, k z ln = 0.1, for 
superlattice period d — 0.251h, d = 0.5/fr, and d^> Ih- 



V. AN INFINITE PERIODIC GRAPHENE 
SUPERLATTICE 

Let us consider an infinite periodic graphene superlat- 
tice consisting of 2D layers parallel to the xy-pl&ne and 
located at z = Id where I — 0, ±1, ±2, ••• , ±oo and d 
is the period. The layers are embedded in a medium 
with background dielectric constant Sb- The dispersion 
equation may be calculated in RPA in the same way 
described in Refs. [HI, [13]. It may be shown that the 
dispersion relation for magnetoplasmons in a superlat- 
tice is obtained by solving 1 — V c (q)U(q,uj)S(q,k z ) = 0, 
where H(q, oS) is the polarization function for a single 
graphene layer defined by Eqs. {!])- ([2]). Also, S(q,k z ) 
is the structure factor determining the phase coherence 
of the collective excitations in different layers given by 
S(q,k z ) = sinh(qd) / (cosh(qd) — cos(fc z d). Note that the 
periodicity ensures that S(q,k z ) is independent of the 
layer index I. Also, the effective-mass model is employed 
to represent the low-frequency electron band structure 
of the layered graphene. At small separations, the low- 
energy bands may be modified by the interlayer atomic 
interactions. In this case, the Landau levels may be dis- 
persive in the k z wave vector. 

We have solved the magnetoplasmon dispersion equa- 
tion in the complex frequency plane. We present only 
the imaginary part of the solution for various values of 
d/ls in Fig. II for k z l H = 0.1. 

The results of our numerical calculations for an infinite 
graphene superlattice show that there are magnetoplas- 
mon modes and modes independent of the wave vector 
corresponding to single-particle transitions between Lan- 
dau levels. Due to a magnetic field, the group velocity 
is negative, as seen over a given range of wave vectors. 
This is analogous to the magnetoplasmons in a single 
graphene layer. Energy transfer from a charged parti- 
cle to the collective modes occurs only when the charged 
particle's velocity has the same value as the phase veloc- 
ity of the collective mode. The most energetic collective 



mode increases for small wave vectors qln < 1 and is 
Landau damped for large wave vectors. It is shown in 
Fig. dj that the imaginary part of the collective mode 
frequencies responsible for the magnetoplasmon instabil- 
ity is appreciably enhanced in an infinite superlattice of 
graphene layers compared to the single layer. Both the 
real and imaginary parts of magnetoplasmon frequencies 
are much larger in a graphene superlattice than in a sin- 
gle graphene layer due to the superposition of the collec- 
tive modes corresponding to oscillations from different 
layers occurring in-phase. The amplification of the col- 
lective mode frequencies increases when a/lu decreases. 
According to Fig. [4] it is clearly shown that the amplifica- 
tion of magnetoplasmons is increased as a/In is reduced. 
When a/In = 0.25 the magnetoplasmon frequencies are 
about twice as large compare to these at a/In —> oo. For 
a/In = 0.5, these corresponding frequencies are larger by 
a factor of 1.5 relative to the result when a/In — > oo. 

While the 2D energy band is not suitable for de- 
scribing the low-frequency electronic properties of bulk 
graphite [33| . the calculated magnetoplasmon frequen- 
cies obtained from our superlattice model are valid. This 
is the case because the separation between neighboring 
graphene layers in the superlattice is much larger than in 
bulk graphite. In a graphene superlattice, the distance 
between graphene layers can be sufficiently large, e.g., as 
assumed in Fig. Ad — 0.25ijj, d — 0.51h, or d — Ih (e.g., 
Iji = 66A at B = 15T), which is large compared to the 
distance between carbon layers in bulk graphite which is 
c/2, where c/a = 2.802 with a — 2.566A denoting the lat- 
tice constant. The significance of this enhanced magne- 
toplasmon instability in superlattices of graphene for de- 
vice applications may lie in possibly utilizing the energy 
of the amplified plasma modes for electromagnetic radia- 
tion in the THz regime, leading to a potential new source 
of radiation based on superlattices of graphene layers. 
For example, with an applied magnetic field B = 10T, 
corresponding to filling factor v = 1, the magnetoplas- 
mon frequency is about 3.6THz. Moreover, the advan- 
tage of such sources of THz radiation is the fact that the 
frequencies corresponding to magnetoplasmon instability 
leading to THz electromagnetic radiation decrease when 
applied magnetic field increases and the parameter a/In 
decreases which results in the possibilities of controlling 
the THz radiation frequencies by changing the applied 
magnetic field. 



VI. DISCUSSION 

We emphasize the appearance of a magnetoplasmon 
instability in a single graphene layer even without the 
application of an in-plane current driving the charge car- 
riers. This instability corresponds to the finite imagi- 
nary part in the frequency of the collective excitations in 
Fig. [TJ There is a plasmon instability in a bilayer semi- 
conductor without an-in-plane current appears only in a 
very small region (compared to the Fermi wave vector) of 
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the wave vector [34J. This difference in the spectrum of 
collective excitations in graphcnc structures compared to 
layered semiconductors is caused by the screening prop- 
erties of the dielectric function in graphene [Icl |35| | and 



2D semiconductors [2~i| . 
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